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A STUDY OF THE EFFECTS OF REYNOLDS 
NUMBER AND MACH NUMBER ON CONSTANT 
PRESSURE COEFFICIENT JUMP FOR SHOCK-INDUCED 
TRAILING-EDGE SEPARATION 

By 


Atlee M. Cunningham, Jr., and Gregory S. Spragle 


SUMMARY 

A study was conducted to address the influence of Mach 
and Reynolds numbers as well as airfoil and planform geometry 
on the phenomenon of constant shock jump pressure coefficient 
for conditions of shock-induced trailing-edge separation 
(SITES) . It was demonstrated that the phenomenon does exist 
for a wide variety of two- and three-dimensional flow cases 
and that the influence of free stream Mach number was not 
significant. The influence of Reynolds number was found to 
be important but was not strong. Airfoil and planform 
geometric characteristics were found to be very important 
where the C jump was shown to vary with the sum of (1) 

r 

airfoil curvature at the upper surface crest and (2) camber 
surface slope at the trailing edge. It was also determined 
that the onset of SITES could be defined as a function of 
airfoil geometric parameters and Mach number normal to the 
leading edge. This onset prediction was shown to predict the 
angle of onset to within +1° accuracy or better for about 90% 
of the cases studied. 
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INTRODUCTION 


The phenomena of shock boundary-layer interaction and 
shock-induced separation as well as their influences on 
aircraft performance have been the subjects of intense 
research of many years. Shockless and other airfoil design 
techniques have minimized the adverse effects of these 
phenomena on transonic cruise vehicles. However, fighter 
aircraft which maneuver at transonic speeds frequently 
encounter extensive shock- induced separations due to the high 
incidence required to achieve high normal force. Optimum 
maneuver capability, which is of upmost importance to fighter 
aircraft survivability, is currently developed on the basis 
of limited experimental studies. Such design techniques are 
expensive and time consuming. More importantly, true optimum 
designs are most likely never realized because of the 
designers limited visibility and configuration inflexibility. 
Therefore, expansion of the design data base with theoretical 
or semi-empirical analysis and design methods is needed in 
order to achieve more optimum designs. 

Development of analytical methods for treating 
shock-induced separation requires a basic understanding of 
the phenomenon. Pearcey provided a very illuminating 
discussion of the process for turbulent boundary layers in 
transonic flow over airfoils in Reference 1. His criteria 
for onset of shock-induced separation stated that it would 
occur when the shock pressure ratio reached 1.4 and the 
downstream pressure reached sonic value. However, he did not 
address the nearly constant pressure-rise amplitude through 
the shock once separation reached the airfoil trailing edge. 
Cunningham (Reference 2) , found that if this pressure rise on 
the airfoil surface was cast as a pressure coefficient jump, 
<C >, a nearly constant value resulted regardless of the 

r ( 

upstream shock Mach number, angle of attack, or span statxon 
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location for finite wings. The value of <C p > was found to 

vary from 0.44 to about 0.67 with Reynolds number for the 

limited cases examined. Thus, it appeared that, depending on 

how energetic the boundary layer was, only a given amount 

of the free-stream dynamic pressure would be recovered across 

a shock terminating on that boundary layer. The boundary 

layer would thicken aft of the shock to form the wedge angle 

necessary to satisfy the oblique shock relations as pointed 

out by Pearcey. The free surface defined by the thickened 

boundary layer was therefore perceived to couple with the 

lower surface flow in a universal manner analogous to the 

Kutta condition for attached flow, such that the limit value 

of <C > was maintained. 

P 

As a result of the above findings, the study summarized 

in this report was conducted to determine the effects of 

Reynolds number and Mach number on the value of <C p > f° r 

shock induced trailing edge separation. The study addressed 

the influence of airfoil geometry on these effects in 

two-dimensional flows as well as planform geometry for 

three-dimensional flows. The first task was to develop an 

appropriate data base from published or other available 

sources of wind tunnel data with Reynolds and Mach number 

variation sufficient to establish constant <C > trends with 

P 

Mach number, Reynolds number, and the geometric parameters. 

The next task involved processing the experimental data into 

the "shifted C plot" format from which the <C > values were 
p P 

determined. Finally, the results were analyzed to establish 

appropriate parameters that might be used by the designer or 

analyst to establish constant <C p > values for a new aircraft 

design as well as conditions under which the phenomenon might 

occur. These results will also be helpful in evaluating wind 

tunnel test results and how they might be altered under full 

scale conditions. 
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SYMBOLS 


A ted 


b 

c 



Non-dimensional parameter for indicating the 
transition to shock induced trailing edge 
separation (Equation 3) 

Wing span, in 

Wing or airfoil chord, in 

Pressure coefficient (p - p^ )/q 

Value of Cp for sonic velocity 



Pressure coefficient value just forward and aft 
of the shock respectively 


‘V 

< s > 


* 


Shock jump Cp value, Cp 2 - C ^ 

Reduced <C p > for airfoil geometry effects 
(Equation 2) 


Transonic similarity parameter modified to 
account for airfoil incidence (Equation 1) 


M 


M 


00 




n 


MAC 

P 


Mach number 

Free stream Mach number 

Mach number normal to the leading edge, 
Moo ( Cos A le ) 

Mean Aerodynamic Chord, in 

Local static pressure, psi 
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Free stream static pressure, psi 

q Free stream dynamic pressure, psi 

c/R Airfoil curvature at the upper surface crest 

non-dimensionalized by local wing chord 

R e Reynolds number based on wing chord for 

two-dimensional flow and MAC for 
three-dimensional flow 

t Wing or airfoil thickness, in 

x Chordwise coordinate, in 

y Spanwise coordinate, in 

( x /c) crest Chordwise location of the upper surface crest 

(l~x/c) cres ^. Chordwise location of the upper surface crest 
relative to the trailing edge 

a Angle of attack, positive nose up, deg. 

Wing twist angle, measured at the 2/3 span 
station, deg.* 

a ted Angle of attack at which the trailing edge 

pressure diverges and switches from positive to 
negative signaling the onset of shock induced 
trailing edge separation, deg. 

* 

Angles are defined in degrees except in equations 1, 2 and 
3 where they are defined in radians. 


err 


a 


a ted - a ted , , deg. 

meas pred' J 


One half airfoil thickness ratio, (t/2c) 


§lef,5tef Deflection of leading and trailing edge flaps, 

degrees (positive nose/trailing edge down) , deg. 


§tec /DELTEC slo P e °f the airfoil camber surface at the 

trailing edge (positive trailing edge down) , 
deg. * 


§ tel ,DELTEL Slope of the lower airfoil surface at the 

trailing edge (positive trailing edge down) , 
deg. * 


Y Ratio of specific heats, 1.4 for air and nitrogen 

A le Leading edge sweep angle, deg. 

17 Span station as a fraction of wing semi-span, 

y/(b/2) 


Angles are defined in degrees except in equations 1, 2 and 
3 where they are defined in radians. 
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BACKGROUND 


The development of shock- induced separation that extends 
to the trailing edge is a continuous process that can be 
separated into distinct phases, as discussed by Pearcey in 
Reference 1. These phases can be described as (1) the 
initial phase, where the shock is weak and the separated 
region is small and confined to the foot of the shock; (2) 
the transition phase, where the separation bubble grows 
rapidly until reaching the trailing edge; and (3) the final 
phase, where the trailing edge remains separated. 

The pressure data shown in Figure l on a supercritical 

airfoil in two-dimensional flow (Reference 3) is a good 

example of typical pressure distributions that occur in the 

vicinity of the shock interacting with a separated turbulent 

boundary layer during the first phase. At =0.71°, 1.18°, 

and 1.567°, it will be noticed that the pressure coefficient 

just aft of the shock where the curves break (x/c-0 . 62-0 . 65) 

is about C * , which is the sonic value of C . Instead of 
. P P 

trying to follow Rankine-Hugoniot conditions across the shock 

for increasing a , the velocity is remaining constant at the 

sonic value. 


Schlieren photographs of typical shock/boundary-layer 
interactions, where the boundary layer is turbulent and 
separates at the shock, show that the shock is not normal to 
the airfoil surface. Comparing the shock angles and the 
"wedge" angles produced by separated boundary layers, the 
numbers are typical of supersonic wedge flows with oblique 
shocks where M = 1.0 is produced behind the shock. For 
example, for an upstream Mach number of 1.126, y = 1.4, and a 
wedge angle of 2°, the shock angle is 18° from the normal and 
the Mach number aft of the shock is 0.9949 (Reference 4) . As 
another example, for a wedge angle of 6° and an upstream Mach 
number of 1.285, the shock is inclined 24° from the normal 
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and the downstream Mach number is 0.9915. 


Considering these observations and the fact that the 
flow is subsonic in the boundary layer, it makes sense that 
the local flow in the vicinity of the shock/wall intersection 
can adjust itself in response to a flow disturbance in order 
to maintain a velocity slightly less than sonic just aft of 
the shock. The maintenance of near sonic velocity is logical 
since it is the highest subsonic velocity that will permit 
forward propagation of disturbances. Thus we are led to the 
conclusion that an appropriate boundary condition to be 
satisfied across a shock on the wing surface is the 
maintenance of constant C (or velocity) aft of the shock if 

tr 

the boundary layer is turbulent and locally separated. This 
condition was incorporated in the transonic perturbation 
method (Reference 2) and was found to provide excellent 
agreement between theory and experiment for weak shocks on 
finite wings. 

The transition phase usually occurs very rapidly 
according to Pearcey for two-dimensional flows and has been 
observed to do the same for three-dimensional flows. On 
finite wings, the spanwise spreading of shock- induced 
trailing-edge separation is the mechanism by which the 
transition seems to occur at a given span station as shown in 
Reference 3 for the ONERA M-6 wing at M = 0.92. 

When the shock-induced separation reaches the trailing 
edge in the final phase, the shock- jump conditions change to 
another form of limiting condition. One well-known signal of 
this occurrence is trailing-edge pressure divergence in which 
the pressure coefficient at the trailing edge drops 
significantly from a positive value to a negative value for a 
small increase in a . This is illustrated in Figures 2 and 3 
for V - 0.90 and Figures 3 and 4 for 17 = 0.65. Note also in 
these distributions that the character of the pressure 
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1 variation aft of the shock has changed. More importantly, 
the pressure rise across the shock, <Cp>, is about constant 
; for a = 4° and 6° at >7 = 0.90 and a = 6° at V = 0.65. 

I 

1 

This characteristic is highlighted in Figure 5 where the 
Cp distributions for several span stations and angles are 
plotted together in such a manner as to show the constant 
<Cp> at the shock. The circled symbols near the origin 
indicate how the origin of each plot was shifted relative to 
the distribution, a = 6°, jj = 0.65, so as to superimpose the 
shock jump distribution. The fact that the origin shift is 
at an approximately constant C for a constant a indicates 
that Cp 1 and Cp 2 are constant along the span in the separated 
region. (The format shown in Figure 5 will be referred to as 
"shifted C plots" throughout this report) . 

hr 

In addition to <C > being constant for various a values 

P 

and t] stations, the C distributions just aft of the shock 
for some distance are also similar. Thus, the shock- induced 
flow-separation mechanism along the span must be the same 
with it beginning at the wing tip and progressing inboard. 
Since this mechanism involves flow separation, it was 
expected to be sensitive to Reynold's number. 

Shown in Figure 6 is a plot similar to that shown in 

Figure 5 but for a C-141 model at various angles of attack. 

Although these data are all at a single span station, the 

angle of attack varies from 1° to 4°. The shock 

characteristics are similar to those shown in Figure 5 except 

for the <C p> magnitude and the Cp distribution aft of the 

shock. The ONERA M-6 wing data in Figure 5 were obtained for 

a Reynold number based on the MAC of 11.7xl0 6 , and <C p > i s 

seen to be about 0.44. The C-141 data were obtained for 

Re = 20xl0 6 and the <C > is about 0.57. Figure 7 shows a 

P 

limited set of data from C-5A flight test at two angles where 
again the characteristics are similar except that the <C p > is 
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about 0.62 corresponding to R g = 80x10 . Thus there seemed 
to be an increase in <C p > with increasing Reynolds number. 
This corresponds to the ability of an increasingly energetic 
boundary layer to withstand an increasing shock pressure jump 
before separation. 


Conditions under which the shock-induced trail ing-edge 
separation occur, are expected to be a function of , a 
and airfoil characteristics. It was proposed in Reference 2, 
that a relationship based on the well-known transonic 
similarity parameter for thickness could be used to account 
for a as well. Evaluation of this equation for determining 
when shock-induced trailing-edge separation covered the outer 
30% of the span led to inconclusive results. However, some 
modification yielded a form that did seem to be more 
universal, as is given in the following expression: 





{l - M* 


y+i -i2/3 

( a + 5 ) 

2 J 


( 1 ) 


where S is one half of the wing thickness ratio and Y = 1.4 
for air. Examples of application of the above equation to 
various cases in the data base of Reference 2 are tabulated 
below: 



t/2c (at T) = 0.70) 


a ted 

K ia 

Convair 880 

0.0424 

0.80 

8 • 7° (-) 

0.507 



0.85 

1 

0 

CO 

CO 

0.499 



0.89 

00 

• 

-J 

0 

1 

0.477 

Onera M-6 

0.0500 

0.88 

6 . 0° (+) 

0.565 

Wing 


0.92 

5 • 0° (-) 

0.552 

F-lll TACT 

0.0357 

0.854 

10 . 2° ( -) 

0.470 



0.901 

8 • 2° (-) 

0.491 
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The (-) or (+) “ted notat i° n indicates whether separation 

had already occurred (-) at some lower angle or was about to 

occur at a slightly higher angle (+) . It appears that 

K ta < 0,55 i ndicates that the outer third of the wing should 

be separated. However, if the outer 15% of the wing is to be 

used as the indicator to trigger application of the trailing- 

edge separation model, a higher value of k! would be 

ta 

necessary. 

This background forms the basis for a more extensive 
investigation of the constant <C p > phenomenon for shock 
induced trailing edge separation (SITES) . The first step 
will be to develop an appropriate data base using existing 
experimental data from tests in which SITES occurred. Two 
and three dimensional data for both conventional and 
supercritical airfoils are desired in order to establish the 
correct trends with Reynolds number, Mach number and 
geometry . 


DATABASE DEVELOPMENT 

The data used for refinement of the characteristics of 
the constant <C p > concept were obtained from existing sources 
for the purpose of this study. Reference 3 contained a 
wealth of relevant two-dimensional as well as other 
three-dimensional data besides the ONERA M-6 wing results. 
Another vast source of information resided in both published 
and unpublished data from NASA Langley and Ames. General 
Dynamics has also published a significant data base as 
Volumes III through VII of Reference 2. Data for the C-141 
and C-5 are available in reports published by Lockheed for 
both NASA and the Air Force. 
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Data Requirements 


The requirements necessary for the data to be relevant 
to this study are quite specific. The increment of incidence 
change must be on the order of 1° but not more than 2°, so 
that experimental perturbations are not "smeared" by 
excessive shock motion. The range of maximum incidence must 
be such that several data points are available after the 
initial appearance of trail ing-edge separation at some point 
on the wing. Reynolds and Mach number variations are desired 
for a given configuration. Most importantly, the spacing of 
pressure orifices must be close enough to accurately identify 
the shock pressure-rise characteristics and the shock 
location. 

In order to isolate airfoil effects from those due to 
planform geometry, two-dimensional data are needed with 
sufficient Mach and Reynolds number variation for a fixed 
airfoil shape. Three-dimensional data are needed to 
establish planform effects and to distinguish essential 
differences between two- and three dimensional flows. 

Finally, other data from which the occurrence of shock 
induced trailing edge separation can be identified can be 
used to further refine the definition of onset conditions. 
These last data sets do not have to meet the above 
requirements for determining <C > values, but only need to 
contain information relevant to trailing edge pressure 
coefficient divergence. 

Two-Dimensional Data 

The two-dimensional data used in this study were 
obtained from References 3 and 6 through 9 and data which is 
as yet unpublished on the NACA 0012 airfoil. This data was 
obtained from the Langley 0.3-meter Transonic Cryogenic 
Tunnel (0.3-m TCT) by C. L. Ladson. Examples of other data 
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sources considered but not used because of their failure to 

meet the above requirements are listed as References 10 

through 13. Table 1 summarizes the data sets used to define 

the constant <C > characteristics for SITES flows over 

P 

various airfoils. Airfoil types, Mach and Reynolds number 
ranges and wind tunnel test facilities are listed for easy 
reference. Unfortunately, most of the data as well as all of 
the high R data were obtained in the NASA 0.3-m TCT so that 
it was difficult to determine if wind tunnel characteristics 
were important. The availability of data pertinent to this 
study was limited since SITES occurs at conditions which are 
significantly off-design and are thus of little interest in 
airfoil design development. 

The MBB-A3 and ONERA D airfoils listed in Table 1 were 
chosen specifically because data are also available for 
three-dimensional configurations which use these airfoils. 
Hence, these data are helpful in separating two and three- 
dimensional flow effects. 

Three-Dimensional Data 

The three-dimensional data used in this study were 
obtained from References 3, 5, 14, 15 and 16. Table 2 
summarizes the data sets used to define the constant <C p > 
characteristics for SITES flows over various configurations. 
As mentioned above, two configurations had airfoil sections 
which were included in the two-dimensional data sets. The 
ONERA M-6 wing used the ONERA D airfoil section and the MBB 
wing/body used the MBB-A3 airfoil section. 

Most of the applicable data available were obtained from 
wind tunnel tests, however, C-141 and C-5A data were also 
available from flight test. This extended the Reynolds 
number range to 8 0x10 6 with intermediate points at 20, 45 and 
58xl0 6 . The F-16 pressure model data, ref. 15, also included 
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various leading edge flap deflections. Aspect ratios varied 
from 3 for the F-16 to about 8 for the C-141 and C-5A. Wing 
sweeps varied from about 27° for the C-141 and C-5A to 40° 
for the F-16. Test facilities were different for all 
configurations which included five wind tunnels as well as 
flight test. 


Data for Onset Determination 

These data sets were used to further establish what 
conditions were necessary for the occurrence of SITES. Since 
only one or two incidences were needed, one prior to and the 
other after SITES, many of the cases available in the two- 
and three- dimensional data bases which did not qualify for 
determining <Cp>, would qualify for onset determination. In 
addition, the pressure data for the F-lll TACT model given in 
Reference 2 (Vol . IV) , provided onset information for a 
different configuration with a supercritical airfoil section. 


DATA PROCESSING 

Data processing was accomplished in four steps: (1) 
reviewing existing plots and/or tabulated data sets to select 
appropriate data for analysis; (2) replotting all data in a 
consistent format; (3) producing the shifted C plots; and 
(4) determining <C p> values. In order to facilitate the 
investigation, an existing plotting program was modified for 
the purposes of this study. The plot program was coded for 
the General Dynamics/Fort Worth VAX computer system with an 
input format which eliminated redundant input for multiple 
pressure data sets. The capability was also developed to 
permit shifting of pressure data sets to align shock 
location; this was done by hand in prior investigations, see 
ref. 2. Finally, a technique was incorporated in the program 
to calculate <C > values in a consistent manner through the 
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use of least squares fits of the shifted Cp plot data. 

This section will present a description of the data 
processing for producing the shifted C plots and determining 
<C > values using the least squares technique. 

Jr 

Shifted C Plots 
P 

A sample of the final shifted Cp plot format is shown in 

Figure 8. This sample shows three pressure data sets for a = 

1.5°, 2.0° and 3.0° where data for the two higher a sets were 

shifted so as to match the shock jump characteristics with 

the data at a = 1.5°. The plot ordinate and absissa are 

shown for the a = 1.5° data set which is referred to as the 

* 

"REF" data as denoted on the plot. The sonic value, C , is 
also denoted for the "REF" data. Origin shifts for the 
higher a data are denoted as circled symbols for each of the 
data sets. An aft shift means that the shock has moved 
forward and a downward shift means that just upstream of 
the shock, has become more negative. 

A complete set of the shifted Cp plots is available in 
the Appendix. An index of these plots is given in Table 3 
according to configuration, Mach and Reynolds numbers and 
figure numbers. 


Estimation of <C > 

P 


A technique was formulated for automatically calculating 

<C > values in a consistent manner. The technique makes use 
P 

of local least squares curve fits in three user selected 
regions of the shifted C p plots as shown in Figure 9. The 
three regions, upstream, shock and downstream, are determined 
as approximate bounds by the user by specifying the 
appropriate x/c values. Curve fits of up to quartic in x/c 
are used in each of the regions to define mean curves as 
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illustrated in Figure 9. The two intersections with the 

shock fit curve then become the shock C values from which 

P 

<C > is determined. A shock C slope is also defined as 

r r 

<Cp>/Ax where Ax is the "width" of the shock rise as defined 
by the curve intersections. It is not necessary for the 
curve intersections to correspond with the region boundaries; 
however, they were usually quite close to each other. 

A sample plot produced with this technique is shown in 
Figure 10. The plot format is expanded for higher resolution 
to include only the x/c limits specified by the user. Data 
listed in the table in Figure 10 show the <Cp>, <Cp>/Ax, and 
the corresponding curve fit order for each region (1 means 
linear fit, 2 means parabolic, etc.). Curve fits up to 
fourth order could be used in each region. The regions are 
numbered sequentially, 1, 2 and 3 corresponding to upstream, 
shock, and downstream respectively. As shown in this 
example, several combinations of curve fits are used to 
evaluate stability of the fit. In this case, a 3, 2, 3 
(cubic, parabolic, cubic) fit was selected, however later 
analysis of all data showed that 2, 1, 2 (parabolic, linear, 
parabolic) produced more consistent results. In some cases 
where insufficient data points were available for determining 
a reasonable fit, the <C > values were read by hand as done 
in the early investigations. 


RESULTS AND ANALYSIS 

Results for the variation of <C > with Mach and Reynolds 

hr 

number were obtained for both two and three-dimensional 
flows. Mach and Reynolds number effects as well as airfoil 
geometric influences were separated for two-dimensional 
flows. Wing geometry effects were then determined for 
three-dimensional flows which included both wind tunnel and 
flight test data. The onset of shock induced trailing edge 
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separation was determined for all cases considered and was 
cast in the form of a non-dimensional parameter. 

This section will present the two- and three-dimensional 

flow results for <C > variations with various quantities as 

P 

well as results for the onset parameter study. 

Two-Dimensional Constant <C > Characteristics 

P 

The two-dimensional <C p > data generated in this study 
are tabulated in Table 4 and the pertinent airfoil geometry 
data are given in Table 5 . Variation of <C > with Mach 
number is shown in Figure 11. With exception of the data for 
the ONERA "D" airfoil and the NACA 0012 airfoil at M=0.82, 
all data are bounded between a <C p > of approximately 0.5 and 
0.6. These results imply that there is no conclusive trend 
with free stream Mach number. This is not too surprising, 
however, since the value of <C p > seems to be relatively 
independent of the upstream Mach number for both two- and 
three-dimensional flows. 

The variation of <C p > with free-stream Reynolds number 
(based on the airfoil chord) is shown in Figure 12 . These 
results show a trend of increasing value with R e , however, 
the data spread makes it difficult to determine the exact 
trend. The spread in data at R g = 15 and 30 x 10 6 for the 
NACA 0012 and SC (3) airfoils was used to analyze the effects 
of airfoil parameters since these airfoils were significantly 
different. The results of this study are shown in Figure 13 
where the parameter 


<c p > 


<C P > 


C C / R + S tec ] 


2/3 


( 2 ) 
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is plotted against R g . The basis for this parameter lies in 
the non-dimensional transonic similarity parameter, Eqn. 1, 
which was discussed earlier. Inclusion of Mach number in the 
parameter was not successful, as might be expected based on 
the trends shown in Figure 11. The most remarkable results 

shown in Figure 13 are the reduction of scatter at R =15 

& ^ 
and 30 x 10 and the lining up of these points with the R-4 

g 

airfoil data at R = 40 x 10 . 

e 

In order to reduce scatter at constant Reynolds number, 
the values of <C p> at different Mach numbers but constant 
Reynolds number were averaged together. Replotting this 
"averaged” two-dimensional data in the same formats as in 
Figures 12 and 13, yields more clearly defined trends as 
shown in Figures 14 and 15. Figure 14 allows the separation 
of airfoil effects where as Figure 15 emphasises the 
collapsing of the data to a more orderly variation with R e< 

A possible trend which could be suggested from Figure 15 is a 
reduction of data scatter with increasing R e , as indicated by 
the dashed envelope. The limit to which the data appear to 
converge at higher R fi seems to be constant at about 0.88 
which also lines up with several points at lower values of R g 
that include NACA 0012, SC (3) and CAST 10 airfoil data. 

Figures 14 and 15 also show some added data points for 
two different symmetrical bi-convex airfoils at zero 
angle-of-attack. These data were obtained from Reference 17 
for a 12% thick bi-convex airfoil and from Reference 18 for 
an 18% thick bi-convex airfoil. In both cases, the 
shock-induced separation was reached through an increase in 
Mach number instead of an increase in a . The <C p > values 
for the bi-convex airfoils in Figure 14 agree very well with 
each other, however, the <C > values in Figure 15 do not. 
Since the bi-convex airfoils at zero incidence have symmetric 
separated flows, it is quite probable that they cannot be 
classified in the same way as the non-symmetrical cases. 
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This, however, raises a fundamental question as to how would 
the flow transition from symmetric to unsymmetric when SITES 
already exists on a symmetric airfoil at zero incidence and 
the incidence is increased. According to the observations, 
incidence is supposed to have little influence on the value 

of <C > once SITES has occurred. 

P 

Another point that needs to be made about the data shown 
in Figures 11 through 15 is that many of the tests were 
conducted in the NASA Langley 0.3-m TCT. As shown in Table 
1, data for the NACA 0012, SC(3), CAST 10 and DFVLR R-4 
airfoils were obtained in that facility. Noting the spread 
of these data as well as those from other facilities, 
however, it appears that the data source has no particular 
influence on the grouping of results. 

The reason for the <C > scatter at lower R values and 

P e 

the reduction in scatter at higher R fi is attributed to the 
influence of the location of transition from laminar to 
turbulent boundary layer. With larger distances between the 
stagnation point and the transition point at lower R values, 
the effects of tunnel turbulence and other sources of 
unsteadiness would have a greater impact on where the 
boundary layer would actually transition. Thus, a more aft 
transition point would seem to lead to a greater uncertainty 
in the boundary layer development. This hypothesis is 
somewhat justified on the basis of lower scatter in the data 
obtained for the SC (3) airfoil with fixed transition. These 
data are denoted as the solid square symbols in Figure 15. 

In summary, the two-dimensional results show that the 
average value of <C > for SITES tends to increase as the sum 
of (1) airfoil curvature at the upper surface crest and (2) 
trailing edge camber line slope. In this relationship, the 
airfoil curvature is expressed as the reciprocal of the 
radius of curvature non-dimensionalized by the wing chord. 
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Data scatter is present at low Reynolds number values but 

essentially disappears by about R=30xl0 based on the wing 

chord. Above R = 30xl0 S , there seems to be little effect 

® 6 

due to Reynolds number up to 40x10 . No data were available 

6 

above R = 40x10 , hence the trend could not be verified 
e 

beyond this point. A question was also raised on how 
symmetric SITES may differ from unsymmetric SITES and how 
does the flow transition from symmetric to unsymmetric. 


Three-Dimensional Constant <C > Characteristics 

P 


The <C p > data generated in this study for all 
three-dimensional flow cases are tabulated in Table 6 along 
with pertinent configuration geometry data in Table 7. 
Variation of <C > with Reynolds number (based on wing MAC) is 
shown in Figure 16. Because flight test data for the C-141 
and C-5A aircraft are included, the Reynolds number range was 


extended to 80x10 


The same data are replotted as <C p > in 


Figure 17 which also indicate a trend toward less scatter for 
* 6 

R e >3 0xl0 . Collapsing of the <C p > data in Figure 16 to the 

<Cp>* data in Figure 17 is quite convincing evidence of the 

airfoil parameter influence. An average of the C-141 and 

C-5A <C >* data for R >20xl0 S was found to be 1.013, shown in 
P ® 

Figure 17 as a dashed line, where all of these data points 

fall very close to the line with exception of the C-141 point 

at R^ = 58xl0 6 . 
e 


The airfoil parameters given in Table 7 were determined 
at about 2/3 span of the wing. Although some of the wings 
were twisted, only the airfoil parameters were used. When 
the data for two-dimensional flow in Figure 15 are compared 
with the three-dimensional data in Figure 17 , the values of 
<C p >* appear to be in the same range. Because of the 
differences in two- and three-dimensional flows, it was 
initially believed that the two would be incompatible. 
However, combining these data produces the plot shown in 
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Figure 18 in which it is apparent that the data are 

compatible. All three-dimensional data for R <15xl0 6 fit 

e 

within the scatter boundary developed for two-dimensional 

flow. Data for the C-141 and C-5A at higher values of R all 

© 

tend to fall along a constant line as shown in Figure 17. 

This characteristic is similar to the two-dimensional 

characteristic that suggested <C p >* would be independent of 

R_ for R >30xl0 6 . 
e e 

It is apparent in Figure 18 that a conflict exists with 

regard to geometry and three-dimensional effects. The data 
6 

for R e <15xl0 are limited to the ONERA M-6 wing, F-16 
1/9-scale model and the MBB wing/body configurations in 
addition to two values for a C-5A model. Since these data 
fall within the scatter boundary defined by two-dimensional 
flow, the source of the conflict is the fact that the 
"constant" trends defined at higher R @ for both two- and 
three-dimensional flow do not agree. The constant trend for 
three-dimensional flow is mostly defined by C-141 flight test 
data with a single point at R = 80xl0 6 from C-5A flight 
test. 


The determination of geometric parameters for the C-141 

that were used to calculate <C >* is believed to be a 
. P 

potential source of this disagreement. The values of c/R and 
8 tec were calculated based on interpretation of the airfoil 
description given in Reference 19 rather than actual airfoil 
ordinates as was possible for the other configurations. For 
example, if the sum of c/R and S tec estimate for the C-141 
was increased by 10%, the constant value, <C >*=1.013, would 
be lowered to about 0.95. A geometric error for the C-141 
would only affect the absolute value of the data trend and 
would not change the trend itself since all values would be 
divided by the same constant. 
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* 6 

The C-5A <C > value at R = 80x10 is not subject to 
P e 

the uncertainty associated with the C-141 data; thus, a weak 

■ it 

increase of <C > with R^ is indicated by the results shown 
p e 

in Figure 18. Earlier investigations of this phenomenon 

discussed in Reference 2 considered the possible increase of 

<C > with R 1 / 5 which was based on the variation of a 

P e * 

turbulent boundary layer thickness. Hence, dividing <C > by 

1/5 . P 

R ' should produce a constant trend at high R e « A plot of 

this function for the data in Figure 18 is shown in Figure 

19. The C-141 data do not show a constant trend, however, 

6 6 
the C-5A value at 80x10 and the R-4 value at 40x10 are 

nearly equal. Again, for R >20x10 , the C-141 data seem to 

be displaced which if lowered would tend to fall in line with 

the other high R^ data. The rise of <C >*/R 1//5 for 

6 e g e 

R <30x10 is due to the trend that <C > is centered about a 

e p 1/5 

constant value of 0.88 which, when divided by R ' tends to 

e 

infinity as R g goes to zero. 

Onset of Shock Induced Trailing Edge Separation 

The onset of conditions under which the constant <C > 

P 

phenomenon exists also coincides with the divergence of 
trailing edge pressures. This trailing edge divergence, 
where the pressure coefficient becomes negative signaling 
loss of pressure recovery, is a well know indicator of buffet 
onset. Identification of the onset is of upmost importance 
for design information and represents a significant part of 
this research program. 


The results of the onset study are shown in Figure 20 
where the parameter 


M n l 1 ‘ M n [( 1 >crest ( c/R + S tel >] 


1/3 


ted 


[t( 


c/R + 


° ted + °t + 0,8 S tec 


2/3 


(3) 


( 1 ” X//c ^ crest 

is plotted as a function of Reynolds number for all two- and 
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three-dimensional data. This quantity is also based on the 
transonic similarity term, with several modifications 

which were determined through an extensive parameter study. 
For two-dimensional flow, the wing twist angle, a. , is 
naturally zero but in three-dimensional flow, it is taken at 
the 2/3 span station where other airfoil geometric data are 
determined. The onset study results are tabulated in Tables 
8 and 9 for two- and three-dimensional data respectively. 

The data shown in Figure 20 represent three 
possibilities: (1) before SITES (open symbols), (2) at SITES 
(half-filled symbols) and (3) after SITES (solid symbols) . 
Three-dimensional data are denoted by flagged symbols at 
which onset is considered to occur when the outer 1/3 of the 
wing is in SITES flow. These results show that the 
transition to SITES is indicated approximately by a mean 
value, = 0.19, where below that value the constant <Cp> 

phenomenon should occur. Generally, the spread of transition 
points about the mean value of 0.19 is within the +5% bounds 
(ie. 0.18 to 0.20) with only two points lying outside these 
bounds. There does not appear to be a clearly defined trend 
with Reynolds number except that the average A ted might be 
slightly higher at about 0.195 as Re-*>0. 

The use of Mach number normal to the leading edge for 
three-dimensional data, 

M n = M * cos A ie 

was found to provide the best correlation of results. The 
average of leading edge and trailing edge sweep angles was 
also considered, however, since most of the wings had nearly 
the same sweep for both leading and trailing edges, this 
variation was not significant. The F-16, with A = 40° and 
A te = 0° did show a considerable sensitivity to this effect 
and hence was the deciding factor in selecting the leading 
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edge sweep for determining M, 


n 


The validity of using A^ e( j = 0.19 as an indicator for 
signaling transition to SITES was evaluated by comparing a te( ^ 
values calculated from A ted = 0.19 with actual measured 
values. These data are also listed in Tables 8 and 9 and 
plotted in Figure 21. The solid line shown in Figure 21 
represents perfect correlation whereas the dashed lines 
represent a difference of +1° between measured and predicted 
a ted va l ues * The distribution of measured transition points 
relative to the +1° error band, which represents a reasonable 
accuracy, is listed below: 


Error Band 


No. of Transition Points 


a >1.0 
err 

-1.0°< a <1.0 


a <-l • 0 
err 


grr 


2 

18 

1 


where 


= ated^^ = ^ - ated„„ . 
err meas pred. 


For conditions prior to transition the distribution is 


Error Band 


No. of Pre-Transition Points 


^ a meas 
^ a meas 


a ted 
a ted 


pred 

pred 


) > 1 . 0 ° 
)< 1 . 0 ° 


1 

12 


and for conditions after transition the distribution is 


Error Band 


No. of Post-Transition Points 


^ a meas 
^ a meas 


a ted 
a ted 


pred 

pred 


) >-1.0 
)<-1.0 
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15 

1 


These results provide an evaluation of the accuracy of the 
prediction method on a statistical basis assuming an error of 
+1.0°, where the predictions for transition were accurate to 
within +1° for 18 out of 21 samples. Flow had not 
transitioned for angles less than ( a ted pred +i°) for 12 out 
of 13 samples and it had already transitioned by 
( ateci p re( j“ l0 ) for 15 out of 16 samples. 

The onset study has thus provided a parameter for 
predicting the onset of SITES for a fairly wide range of 
configurations and conditions. The accuracy of the predicted 
angle for transition is nominally +1° although the accuracy 
is probably much better and nearer to +0.5°. The parameter, 
A ted' kased on ly on airfoil and planform geometry 
information and the free-stream Mach number, thus for the 
cases studied, no information was needed from pressure 
distributions . 


Influence of Control Surfaces 


Since the constant <C > phenomenon has been so closely 

Jr 

tied to airfoil geometry in this study, one would expect it 

to be significantly influenced by the deflection of leading 

and/or trailing edge control surfaces. Considering the 

* 

geometric parameters in either <C p> cr A , it would appeal 
that a leading edge flap would not have any significant 
effect. This conclusion was verified with the F-16 data set 
where test conditions included leading edge flap settings of 
0°, 5° and 10°. Data are shown in Figure 22 for S lef 's of 
0°, 5° and 10° at M=0.9, in shifted C p plot format. As can 
be seen, the deflected flap data fall right in with the 
undeflected data. This same result was obtained for the 
other Mach numbers and in no case could an influence on 
either <C p > amplitude or onset be defined. 
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With regard to trailing edge flap deflection, 8 


tef' 


no 


data were available in the data base surveyed. Considering 
the importance of 8 and 8 +^ in the <C_>* and A, 3 , this 


tec 

is unfortunate since 8 


tef 


tel 

would be added to both 


ted' 

8 tec and 


5 tel - For a constant <C p > / trailing edge flap deflected 

down would lead to a higher value of <C >. For a constant 

P 

A ted' the same downward flap deflection would lead to a lower 

value for a ted * Remembering that S tel is in the numerator 

°f A ted , the influence on a ted would not quite be as great as 

on <C > . 

P 


The effect of leading and trailing edge control surface 
deflection would be another area for further investigation. 
The influence of Reynolds number could be more important 
where regions of high curvature exist in the vicinity of the 
hinge line. The control surface effect would also be 
influenced by the control deflection angle as well as both 
two- and three-dimensional flows. For the latter, the 
stabilizing effect of spanwise flow would be expected to be 
very important. Therefore a test program could be envisioned 
that included several airfoils plus leading and trailing edge 
flaps as well as a finite wing planform using at least one of 
these airfoils, all of which would be tested over a wide 
range of Reynolds number. 


CONCLUSIONS AND RECOMMENDATIONS 


A study was conducted to address the influence of Mach 
and Reynolds numbers as well as airfoil and planform geometry 
on the phenomenon of constant shock jump C , <C >, for 
conditions of shock induced trailing edge separated (SITES) 
flows. It was demonstrated that the phenomenon does exist in 
a wide variety of two- and three-dimensional flow cases and 
that the influence of free stream Mach number was not 
significant. The influence of Reynolds number was found to 
be important in that for R e <30xl0 there was more scatter in 
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the <C > results and for R >30x10 there seemed to be a trend 
P e 

for <C p > to be weakly dependent on R e> 


The effect of airfoil and planform geometry was found to 
be very important. It was shown that <C p> increased as the 
sum of (1) airfoil curvature at the upper surface crest and 
(2) camber line slope at the trailing edge was increased. It 
was also determined that the onset of shock induced trailing 
edge separation could be defined in both two- and three- 
dimensional flows as a function of airfoil geometric 
parameters and Mach number normal to the leading edge. The 
resulting prediction technique provided estimates of the 
onset within an error bounds of +1° or better for about 90% 
of the cases studied. 


As a result of the findings of this study, several 
recommendations are offered to answer some of the questions 
that have been raised. The two-dimensional characteristics 
have been shown to be strongly related to the three- 
dimensional characteristics, this enables airfoil tests to be 
conducted to provide a better understanding of 
three-dimensional flows. The recommendations are as follows: 


1. It is recommended that a symmetric airfoil such as 
the NACA 0012 be re-tested starting with a= 0° but with Mach 
numbers high enough to produce SITES. Once SITES is 
established, the incidence should be increased up to the 

This 

test has the objective of answering the question of how does 
a SITES flow transition from a symmetric to a non-symmetric 
flow and what happens to <C > in the process. 


point where the constant <C > phenomenon disappears. 

P 


2. A series of tests is recommended where all or some 
of the four basic airfoils, NACA 0012, SC (3) , CAST 10 and 



DFVLR R-4 , be re-tested to a much higher Reynolds number, 

g 

hopefully 80x10 . This series of tests will have the 

» • • • • it 

objective of determining if the trend of increasing <C > 

. 6 P 

with R^ for R >30x10 holds true for two-dimensional flows, 
e e 

3 . A third test is recommended where one of the above 
four airfoils is used on a finite wing with perhaps several 
leading edge sweeps of say 0°, 20° and 30°. This test will 
provide information on spanwise flow effects at higher 
Reynolds numbers than were available in the data base 
developed for the current study. If the NACA 0012 airfoil is 
used, this test could also include the investigation of 
transition from a - 0, recommendation 1. 


4. Finally, a series of test is recommended where 
leading and trailing edge flaps are investigated in 
two-dimensional flows. In addition, this should be extended 
to three-dimensional flow again using the NACA 0012 airfoil 
suggested in the first and third tests. The objective of 
this series of tests would be to identify the effect of flap 
deflection angles and Reynolds number interplay on the <C^>* 
and A ted parameters. 
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TABLE 1 - DATA SETS FOR TWO-DIMENSIONAL SITES FLOWS 
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TABLE 3 


- INDEX OF SHIFTED C PLOTS 

P 


CONFIGURATION 

MACH 

R e 

FIGURE 

NACA 0012 Airfoil 

0.76 

6X10 6 

A-l 

if 

0.80 

II 

A- 2 

ii 

0.82 

II 

A- 3 

if 

0.74 

9X10 6 

A-4 

ii 

0.76 

15X10 6 

A- 5 

ii 

0.78 

II 

A- 6 

ii 

0.80 

II 

A-7 

ii 

0.74 

30X10 6 

A-8 

SC ( 3 ) -07 12 (b) Airfoil 

0.78 

ioxio 6 

A-9 

ii 

0.79 

II 

A-10 

ii 

0.80 

II 

A-ll 

ii 

0.78 

15xlO S 

A-12 

ii 

0.80 

II 

A-13 

ii 

0.78 

30X10 S 

A-14 

ii 

0.79 

II 

A-15 

ii 

0.80 

II 

A- 16 

ii 

* 

00 

• 

o 

ioxio 6 

A-17 

ii 

* 

0.78 

15X10 S 

A-18 

ii 

* 

00 

r- 

# 

o 

30X10 S 

A-19 

CAST 10-2/DOA Airfoil 

0.78 

4X10 S 

A-20 

ii 

0.79 

II 

A-21 

it 

0.80 

II 

A-22 

ii 

0.80 

6X10 S 

A-23 

ii 

0.79 

IOXIO 6 

A-24 

ii 

0.80 

II 

A-25 

ii 

0.78 

15X10 S 

A-26 

it 

0.79 

II 

A-27 

ii 

0.80 

II 

A-28 


Fixed Transition 
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TABLE 3 (cont'd) 


CONFIGURATION 

MACH 

R e 

FIGURE 

DFVLR R4 Airfoil 

0.78 

4X10 6 

A-29 

it 

0.78 

15X10 6 

A- 3 0 

ii 

0.78 

40X10 6 

A-31 

CAST 7 Airfoil 

0.76 

6xlO S 

A- 3 2 

MBB-A3 Supercritical 

0.80 

6X10 S 

A-33 

Airfoil 




ONERA D Airfoil 

0.81 

4.5xlO S 

A-34 

II 

0.84 

II 

A-35 

ONERA M-6 Wing 

0.88 

12X10 S 

A-36 

ii 

0.93 

12xl0 S 

A-37 

F-16 1/9-Scale Model 

0.90 

2.5xlO S 

A-38 

ii 

0.91 

II 

A-39 

ii 

0.92 

II 

A-40 

ii 

0.93 

II 

A-41 

MBB Wing/Body 

0.90 

1.34xlO S 

A-42 
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TABLE 4. TABULATED C JUMP RESULTS FOR TWO-DIMENSIONAL AIRFOIL DATA 
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READ - The <C > values are determined by manually reading the shifted C plots. 
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READ - The <C > values are determined by manually reading the shifted C plots 
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TABLE 8. TABULATED ONSET RESULTS FOR TWO-DIMENSIONAL AIRFOIL DATA 
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Key to seperatlon Indicator: -1 * Before shock Induced trailing edge separation 

0 ■ At shock Induced trailing edge seperatlon 

1 * After shock Induced trailing edge seperatlon 
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ONERA-D 0.8100 4.5000 0.0000 6.0000 -1.0 0.1651 0.0000 

(TUNNEL 47/R1-CH) 

0.8500 4.5000 0.0000 4.0000 -1.0 0.1672 0.0000 


TABLE 9. TABULATED ONSET RESULTS FOR THE THREE-DIMENSIONAL CONFIGURATIONS 
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Key to seperatlon Indicator: -1 * Before shock Induced trailing edge separation 

0 * At shock Induced trailing edge seperatlon 

1 = After shock Induced trailing edge seperatlon 
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Figure 2 Upper Surface Pressure Distributions on the ONERA M-6 Wing 
M=0.92, a -2 deg. (Reference 3) 
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Figure 3 

Upper Surface Pressure Distributions on 
M=0.92, a =4 deg. (Reference 3) 

the ONERA M-6 Wing 



Figure 5 Shock Jump Similarity for the ONERA M-6 
Wing Upper Surface Pressures at M=0.92, 
a =4 deg. and 6 deg., 1^=11.7x10®. 

(Reference 3) 
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Lgure 8 


Shifted Cp Plot Format, CAST 10-2/DQA 2 Airfoil, M=0.78, 
R g =4 . 01x106 (Reference 7) 
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Figure 9 Scheme for Determining Shock Jump Cp Values from Shifted 
C p Plots, CAST 10-2/DQA 2 Airfoil, M=0.78, Re=4.01xl0 6 
(Reference 7) 
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Figure 12 Influence of Reynolds Number on <C^> for Raw 
Data in Two-Dimensional Flow 
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Figure 13 Influence of Reynolds Number on < Cp > for Raw 
Data in Two-Dimensional Flow 
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Figure 14 Influence of Reynolds Number on <Cp> for 
Averaged Data in Two-Dimensional Flow 
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Figure 15 Influence of Reynolds Number on<Cp>* for 
Averaged Data in Two-Dimensional Flow 
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Figure 16 Influence of Reynolds Number on <Cp>for Raw 
Data in Three-Dimensional Flow 
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Figure 17 Influence of Reynolds Number on ^Cp)* for Raw 
Data in Three-Dimensional Flow 
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Influence of Reynolds Number on 
<Cp)* for all Two- and Three- 
Dimensional Flow Data 
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Figure 19 Influence of Reynolds Number on 

^p^*/^ 0,2 for all Two- and Three- 
Dimensional Flow Data 
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Figure 20 Influence of Reynolds Number on the SITES Onset 
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Figure A-5 NACA 0012 Airfoil, M=0.76, R e =15xl0 6 
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Figure A-23 CAST 10-2/DOA Airfoil, M=0.80, R 
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Figure A-24 CAST 10-2/DOA Airfoil, M=0.79, R =10x10 
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Figure A-25 CAST 10-2/DOA Airfoil, M=0.80, R =10xl0 6 
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Figure A-32 CAST 7 Airfoil, M=0.76, R e =6xl0 6 


74 


Cp ( r e f ) 


(ref) 0 ALPH=2 .01 
□ ALPH=4 . 47 
* ALPHAS. 0 
0 ALPH=4.74- 



5 0.6 0.7 0.8 0.9 

ref) 


Figure A-33 MBB-A3 Supercritical Airfoil, M=0.80, R 
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Figure A-34 ONERA D Airfoil, M=0.81, R =4.5x10 
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Figure A-40 


F-16 1/9-Scale Model, M=0.92, 


R =2,5xl0 6 
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Figure A-41 F-16 1/9-Scale Model, M=0.93, R =2.5xl0 6 
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Figure A-42 MBB Wing/Body, M=0.90, R e =1.34xl0 6 
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